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Abstract. Green's functions are gauge-dependent quantities. Thus, the manifestation of confine- 
ment in these correlation functions also depends on the gauge. Here we use lattice gauge theory to 
study the gluon and the ghost propagators in a gauge (the so-called A -gauge) interpolating between 
the Landau and the Coulomb gauge. Results are compared to the usual minimal Landau gauge. 
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Green's functions encode completely the non-perturbative properties of a quantum 
field theory. Thus, the infrared (IR) behavior of the gluon propagator D(p) in Yang-Mills 
theory should contain information on confinement. In particular, a violation of reflection 
positivity (at large separation) for the gluon propagator is considered an unambiguous 
signal for the confinement of gluons [1]. This issue has been widely studied in Landau 
gauge. Continuum methods, such as Dyson-Schwinger equations find such a behavior 
for the gluon propagator [ 1]. A clear sign of violation of reflection positivity for the 
gluon is also obtained using lattice numerical simulations in the 3d case |0] while, in 
4d, finite-size effects and technical limitations have not allowed yet a verification of this 
result J^. At the same time, the Landau ghost propagator shows a divergent behavior 
[fTJ, [3|] , stronger than p~ 2 for small momenta p. This strong IR divergence corresponds 
to a long-range interaction in real space, which may be related to quark confinement. 
On the other hand, in order to understand the confinement mechanism, it is important to 
verify how these IR features are modified when considering other gauges |4]. Thus, it is 
important to investigate these correlation functions in a wide class of gauges. Particularly 
important are gauges depending on a parameter, since one can interpolate between 
different gauge conditions. This allows one to investigate changes in the correlations 
functions while varying smoothly the gauge parameter (i.e. the gauge condition). 

Let us recall that the gauge-fixed Lagrangian density in the (Euclidean) continuum 
can be written as 

56 = 4?*W?v + + c a d',Dfc h . (1) 

Here, dL = Xdo + dj, is the gauge field, Dff is the usual adjoint covariant derivative, 
F£ v is the field-strength tensor and the ghost fields c a and c a have been introduced in 

order to exponentiate the Jacobian, i.e. the Faddeev-Popov operator d'^D a ^, introduced 
by fixing the gauge. Note that, for A = 1, the quadratic term in dLA" corresponds to 
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FIGURE 1. The gluon propagator D tr (po, \p\) as a function of A for the lattice volume 40 3 at [5 = 4.2. 
In the left panel we report data forD tr (po,0) while in the right panel we consider D tr (0, \p\). We represent 
A = 1 with black squares, A = 1/2 with open circles, A = 1/10 with diamonds and A = 1 /20 with stars. 



the usual perturbative definition of linear covariant gauges with gauge parameter £ . In 
particular, £ = 1 is the Feynman gauge and the limit £ = corresponds to Landau gauge. 
On the other hand, for £ = 0, the parameter A allows to interpolate between the Landau 
gauge (A = 1) and a Coulomb-like gauge (A =0). Note that the usual Coulomb gauge 
d{Af = has to be supplemented by further constraints q] to yield the limit A = 0. 

Here, we present results for the case t, = (the so-called A -gauge) obtained using 
lattice numerical simulations for the SU (2) gauge group. To be able to reach the IR 
asymptotic regime, we consider the 3d case, extending earlier studies of this subject @]. 
Clearly, for A < 1, Euclidean invariance is explicitly broken by the gauge condition. 
In this case, there are two independent kinematic variables, namely po and \p\. At the 
same time, the gluon propagator D^ v (po, \p\ ) also has two independent tensor structures: 
its temporal component Dqo(po, \p\) and its spatial projection D ti (po 1 \p\). Predictions 
exist [5] that both components of the gluon propagator vanish at zero total momentum, 
P 2 o+\p\ 2 = 0. 

In Figured we report results for D ix (po, \p\) as a separate function of po and of \p\ 
for different values of A. We find that, as a function of \p\, the maximum of D tr moves 
to larger momenta as A decreases. On the contrary, when considering its dependence on 
po, we see that the maximum of D tr moves to smaller momenta and that, already for 
A = 1/10, the maximum can no longer be resolved. In both cases, for small momenta, 
the propagator decreases as A gets smaller. The function Z)oo(0 7 \P\) lS considered in the 
left panel of FigureEJ In this case, as A decreases, the maximum of the propagator moves 
to smaller momenta while the propagator increases in the IR. This is somewhat expected, 
since in Coulomb gauge this tensor component diverges O]. Thus, as A decreases, 
the gluon propagator exhibits more and more a Coulomb-like behavior. Finally, the A 
dependence for the ghost propagator is shown in the right panel of Figure |2| We see 
that, as a function of \p\, it is essentially independent of A. On the other hand, when 
Dg(Po, \P\) is considered as a function of po we find that it increases approximately 
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FIGURE 2. In the left panel we show results for the gluon propagator A>o(0, \p\), while in the right 
panel we report data for the ghost propagator D<j(0, \p\). In both cases we consider various values of A 
(symbols and lattice as in Figure^. Note that the quantity Doo {po , 0) vanishes identically due to the gauge 
condition. 



linearly with X [8]. These results support the findings of recent DSE studies Jsfa. A more 
extensive numerical study for this interpolating gauge will be presented elsewhere 

Finally, let us note that the study of the linear covariant gauges on the lattice is 
much more involved. In particular, at £ > and A = 1, the continuum gauge condition 
BA'aipc) = A a (x) cannot be implemented for any finite /3 because A a (x) has a Gaussian 
distribution, while the values of dA^(x) are clearly bounded. Results for the linear 
covariant gauges will be presented elsewhere iJ. 
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